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1. Introduction 



Coherent control of nanoscale devices is one of the main topics of current research on 
electron transport in systems on the nanometer scale. A manifest realization of control 
is given by pump-probe schemes known from molecular physics [l}[2]- In the context 
of electron transport the "pump" and "probe" steps consist in switching to and from a 
regime, where transport is largely blocked (3J. This provides an effective decoupling of 
the electronic system from the reservoirs and allows for a coherent evolution between 
pump and probe triggers. Such a setup has been successfully used to coherently control 
charge [i] and spin qubits [5j|6j in double quantum dots (DQDs). 

The theoretical description of these experiments is a very demanding task, since 
a fully time-resolved calculations are necessary. Although several formalisms exist 
to address this issue, only very few numerical schemes are available in this context. 
Typically, the numerical approaches to time-dependent electron transport with arbitrary 
driving rely either on equations of motion for non-equilibrium Green functions (NEGF) 
[7-12 or generalized quantum master equations (QME) for the reduced density 
matrix 13 , 14 . Such calculations are very helpful to gain a deeper understanding of 
the mechanisms relevant for coherent control, since time-resolved quantities, such as 
occupations and currents, are readily accessible. 

In this article we concentrate on the experiment of Fujisawa and coworkers on 
coherent manipulation of charge states in double quantum dots jij. One of the hallmarks 



of the experiment is the observation of oscillations of the so-called number of pulse- 
induced tunneling electrons as a function of pulse length (see also Sec. 1.1). These 
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oscillations are commonly associated with coherent tunneling processes between the two 
quantum dots during the manipulation stage. While this picture qualitatively explains 
the main features of the experiments, it neglects the influence of the initialization and 
the measurement stages on the coherent time-evolution. As we will show, consistently 
considering the whole pump-probe sequence provides even stronger evidence for actual 
coherent control. Moreover, the new insights may help to gain a deeper understanding 
of the dynamics in nanoscale devices. 

First of all we briefly introduce the concept of charge qubits in double quantum dots 
and review the main results of the coherent manipulation experiment. The theoretical 
model and the relevant tools used in this work are presented in Sec.[2j Numerical 
results and a detailed analysis using an analytically solvable model are given in Sec.[3j 
The article concludes with a summary in Sec.|4} 



1.1. Charge Qubits in Double Quantum Dots 

From quantum information theory it is well known that in principle any two-level system 



may be used as a single bit of (quantum) information, which is then called qubit 15 
In laterally coupled quantum dots one may find such two-level systems by considering 
charge states denoted by (Ne,N T ), i.e., systems with Ni electrons in the left and N T 
electrons in the right dot. An electron tunneling from left to right corresponds to the 
sequence (N t , N r ) ->• (N t +1,N T ) -> (N e ,N r +l) -> (N e ,N r ). The probability for such a 
tunneling event is largest if the two charge states, \£) = (N(+1,N T ) and |r) = (Nt, N T +1) 



have a vanishing energy difference e 16 17 . In the vicinity of such a resonance the 



coupled quantum dots can be described by a two-level system, which is characterized 
by the energy difference e and the interdot tunnel coupling T c |4 , 17 . Correspondingly, 
the associated qubit is called charge qubit. 

The Hamiltonian of the two coupled charge-states in the DQD, i. e. two coupled 
energy levels, reads 

#dqd = £e(t) c\ci + e T (t) c\c v + T c (c\c r + cJq) + Uc\cic\c r . (1) 

The operators (c n ) create (annihilate) an electron in the left (n = £) or right (n = r) 
dot, respectively. The interdot charging energy U suppresses double occupancy of the 
DQD. The time-dependence of the energies and e T may be given by external gate- 
voltages. Figure [T] shows the corresponding energy scheme for a charge qubit additionally 
coupled to a source and a drain reservoir. In order to use a quantum system as a qubit, 
it has to fulfill at least the following three conditions fl8] : (i) initialization of the qubit 
into a well defined state, (ii) application of unitary operations (quantum gates) and (iii) 
readout of the qubit state. Using a pump-probe scheme with rectangular bias-voltage 
pulses, these three steps have successfully been implemented for a charge qubit in a 
double quantum dot El. The scheme reminds in many ways of the usual pump-probe 



experiments with atoms or molecules 19 21 . In the present case the raising edge of 



the pulse triggers the coherent dynamics (pump), while the trailing edge starts the 
measurement (probe). 
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Figure 1. Pump-probe scheme for the coherent manipulation of a charge qubit [4j. 
a) Initialisation (Vsd — V)), b) manipulation (Vsd = V±) and c) measurement phase 
(Vsd = Vq). d) Time dependence of the bias-voltage pulse Vsd(*)- e ) Current (orange 
line) as a function of time. The stationary value Jo is taken before the manipulation 
phase, t < 0, and for long times, t — > oo. The number of transferred electrons N 
(yellow-shaded area), without the stationary portion (gray-shaded), has contributions 
from the manipulation and the measurement phase. 



The main idea of the experiment consists in suddenly switching between the 
Coulomb blockade regime and a transport regime by using a bias-voltage pulse. In the 
former case the DQD is effectively isolated from the reservoirs since sequential tunneling 



is strongly suppressed 22 . This provides the possibility to coherently control the charge 
qubit |3 1 . The transport regime is used to initialize the system and to readout the charge 
state after manipulation. The whole sequence is schematically shown in Fig.[TJ Before 
and after the pulse, i. e. during the initialization and measurement phase, the source- 
drain voltage is Vsd = Vo and transport through the dot is possible since both charge 
states are within the transport window defined by the source-drain voltage. During 
the pulse, i. e. during the manipulation phase, the source-drain voltage is switched to 
Vsd = V\ and no transport is possible. 

In the experiment, the time-averaged current is measured as a function of pulse 
length r and energy difference e = e T — Eg. The latter is tuned by applying a suitable 
gate- voltage V g to the right quantum dot (not shown in Fig.[T|. Notice, that due 
to capacitive couplings [I], the energy difference e for a given gate- voltage V g also 
changes with the source-drain voltage, which is why we use Eq and E\ for high and 
low source-drain voltage, respectively, in Figs.[pL-c. The pulse is repeated with a 
repetition rate f Tep = 100 MHz and by using a lock-in technique the pulse-induced 
current J p is obtained, which does not contain the asymptotic (stationary) current of 
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the initialization and measurement phase. With J p one finally gets the number of pulse- 
induced tunneling electrons N = J p /ef iep [4J. It is then argued that this quantity is 
equivalent to the occupation of the right dot at the end of the pulse |4l. As we will 



show in Sec. 3J2 this assumption is not always valid. The experimental results show 
pronounced oscillations of N as a function of pulse length. Hereby, the frequency and 
the amplitude of the oscillations depend on the energy difference between the charge 
states. These oscillations are interpreted as a signature of coherent tunneling between 
the charge states (Rabi oscillations). Apart from the oscillations there are two other 
noteworthy features in the experimental results: for a fixed pulse length the function N 
is asymmetric around S\ — and, in particular for e = 0, the number of pulse-induced 
electrons takes negative values. It has been argued that these features are artifacts of 
incomplete initialization or imperfect pulse shapes |4j. In the remaining part of the 
article we will demonstrate, that this must not be the case and the additional features 
instead strengthen the view of a coherent manipulation. 

2. Theoretical Description 

After specifying the Hamiltonian we briefly present the time-dependent description of 
the DQD using a quantum master equation. Furthermore we discuss a model which 
allows for deriving analytical expressions and facilitates a detailed analysis. Specifically, 
we consider the large-bias limit and use a description of the dynamics, which has 



originally been developed for photon-assisted transport through DQDs 23 



2.1. Setup 

Apart from the DQD Hamiltonian [Eq. Q], the total Hamiltonian also contains terms 
describing the electron reservoirs and the tunnel coupling, 

H = H D q D + H TCS + H tun . (2) 

The reservoirs are described as usual by non-interacting electrons, 

-f^res — 2^ ^^ £ ak{t)bakbak , (3) 
ct=L,R k 

where {b\ k } and {b a k} are the electron creation and annihilation operators for the a- 
reservoir state k, respectively. The reservoir single-particle energies have the general 
form e a k{t) = e° ak + A a (t) with the A a accounting for a time-dependent bias. The 
tunneling Hamiltonian for the linear setup of two quantum dots in series, which are 
each coupled to a single reservoir, reads explicitly 

#tun = [ T Lk/bi k c e + T KktV b ] Rk c r + h.c. (4) 

k 

= E E l B ™ s n +) + s t ] BL ] ] , (5) 

a=L,R n=£,r 
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with T ak) n denoting the coupling matrix element between QD n = £, r and the fc-th mode 
of the respective reservoir a = L, R. In the second line of the tunneling Hamiltonian 
we have introduced the abbreviations B a V = Y^k^ak,nb ak , £w = Y^k^ak^ak an d 
= c n , = c]j, which will simplify the notation later on. 

Further, making the wide-band assumption renders the tunnel-coupling elements 
independent of the reservoir state, T ak)n = T a>n , which yields a flat spectral density, 

T amn = 2nT a , n T* m ^ 5(e - e° ak ) . (6) 
k 

2.2. Time-local Quantum Master Equation (TLQME) 

The time-evolution of the density operator g characterizing the total system is 
determined by the Liouville- von Neumann equation, 

ih^g=[H,g}_ , (7) 

where H is given by Eq. (|2]). Since we are only interested in the properties of the 
DQD, we can trace out the reservoir degrees of freedom and thus obtain the reduced 
density operator a = Tr rcs g. In order to get an equation of motion for a, we employ 



the time-convolutionless projection operator technique 24 -26 , which yields a time-local 
differential equation for a [13], 



ij^(t) = [H s ,a(t)]_ -i££ (js£\A£l(t)a(t) - <r(t) A«(t) 



a m 



S^\K^{t)a{t)-a{t)N^{t)] . (8) 



+ 

Here and in the following we set K — 1. The TLQME is supplemented by the following 



auxiliary operators 13 



A&(«)= £ J dt'C^\t,t') Ul(t,t')S^U s (t,t') , (9a) 



' t 

t 



AS(t)= £ Jdec^it^uiit^s^usit,^. (9b) 



to 



These auxiliary operators are determined by the free DQD propagator Us(t,t') and the 
reservoir correlation functions Cam (t,t'). The former is given in terms of the DQD 
Hamiltonian, Ug(t, t') = Texp ^— i f t , dt"H D q D (t") j , where T is the usual time-ordering 
prescription. The correlation functions describe the influence of the reservoir on the 
system dynamics. In the current context they are given by [13] 

C<fc% = / ^ amm f a (e) exp (i jf dt" [e+A^t")]) , (10a) 
^L +) (t, = / Jr amm [l-/ a (e)] exp f-i jf eft" [e+AaCO]) • (10b) 
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The Fermi function f a (e) characterizes the initial state of reservoir a, and the spectral 
density T amm has been defined in Eq. ^ . In order to include the energy-level broadening 
induced by the back-action of the reservoirs on the DQD, we modify the correlation 
functions by multiplying with an exponential decay factor 27^29 . Consequently, we 
have C { a 2\t,t') oc exp(-r amm (*-t')/2) 



2S 



Although the equation of motion for the reduced density matrix is local in time, 
it is still very demanding to solve it numerically An efficient method to propagate a 
is the so-called auxiliary-mode expansion technique, which is based on a decomposition 
of the Fermi function [13]. As a consequence the operators and are also 
decomposed into a sum and for each term an individual equation of motion can be 



derived. The details of this procedure are given in Appendix A The auxiliary-mode 
expansion has successfully been used in the context of time-nonlocal and time-local 



quantum master equations 13 , 14 , 30 and has recently been applied to non-equilibrium 
Green functions fT2l. 



It remains to discuss the calculation for the time-resolved electric current through 

The electric current 



the tunneling barriers in a way consistent with the QME 31 



through the tunneling barrier a is given by the rate of change of the particle number in 
reservoir a 27,32], 

d 



Ja(t) 



e-(N a } = -ie([H,N a ]} 



2elm 



Plugging the formal solution of the Liouville- von Neumann equation [Eq. (m)] in the 



interaction representation into Eq. ( 11 ) and keeping only terms up to 2nd order in H tun , 



one finds within the time- local approximation for the time-dependent current |13| , 

J a (t) = 2eRe ^Tr DQD {c m AW(i)a(i) - c m a(t)KW(t)} • (12) 



This expression is consistent with the TLQME, which is also of second order in H t , 
Note, that in order to calculate the current one only needs the auxiliary operators A 
Aim and the reduced density operator a. 



(+) 

am , 



2.3. Markovian Approximation in the Large-Bias Limit (LBL) 

Considering the experimental situation, one is lead to an even simpler description of 
the dynamics in the DQD. Firstly, for a large inter-dot interaction strength U only 
the following three states are relevant: |0) , \£) = c\ |0) and |r) = c\ |0) 



23 



These 



correspond to an empty DQD, one electron occupying the left dot and one electron 
occupying the right dot, respectively. These three states are used as a basis for the 
reduced density operator in the following considerations. The second simplification 
arises due to the large bias in the experiment. In this case, using the Markovian limit 



of the QME [Eq. ®] is well justified 23,33 . 
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For the initialization and measurement phase (cf. Fig.JT]) one obtains the following 



differential equations for the components of the reduced density matrix 23 , 34 



&ee = — iT c (cr&— a r e) + T L (1— a«- a„) , (13a) 
<x rr = - iT c (a T £-(j£ T ) - r R cr rr , (13b) 
a h = - iT c (a u -a TI ) + ie a £l - (T K /2)a b . (13c) 

Obviously, the coupling to the reservoirs introduces transitions between the charge 
states |0), \£) and |r). The first two equations describe the evolution of the charge- 
state occupations, which change due to tunneling between the two dots and because of 
tunneling from and to the reservoirs. The remaining equation yields the dynamics of 
the coherences a r £ = a* b . Tunneling out of the DQD leads to loss of coherence. 

In the manipulation phase an effective decoupling of the DQD from the reservoirs 
is achieved by switching the source-drain voltage. An electron can still enter the DQD, 
but leaving the system is strongly suppressed. In contrast to the initialization phase, 
where tunneling out of the DQD leads to a vanishing coherence, here the dynamics stays 
approximately coherent. Therefore, we assume the following equations of motion during 
the voltage pulse, 

& u = - iT c (<T b -(T t i) - 7<r« + T L {l-a u -a n ) , (14a) 
cr rr = - iT c {a H -a h ) - 7<x rr + T R (l-a ei -a„) , (14b) 
&t T = - iT c (a«-cr rr ) + i£icx & - ja iT . (14c) 

The rate 7 has been introduced to account for additional processes leading to 
decoherence, such as background-charge fluctuations or back-action of the reservoirs 
on the DQD. In principle, these processes are also present in the other stages, but there 
they are dominated by the transport from source to drain. The expressions for the 
time-resolved currents are very simple and may, for instance, be extracted from Eqs. 
(13) and (14). They are explicitly given by 

J R (t) = er R cr rr (t) , (15a) 
J R (t) = -eT R [l-a u {t)-a rt {t)} +e 1 a rr {t) , (15b) 



in the initialization/measurement and manipulation phase, respectively 23 33 



3. Results and Discussion 



In the following we will discuss a DQD system with parameters based on the 
experimental values |4j. We summarize all relevant quantities in Table 1. As sketched 
in Fig.[TJ we consider a perfect rectangular pulse with duration r and assume that the 
DQD is an stationary state at t — 0, when the manipulation phase starts. In the 
numerical calculation, at t — the source-drain voltage switches from Vq to V\ and at 
t — r it switches back. The level of the left dot is fixed at Ei = — eVo/2. In order to get 
an energy-resolved picture we can shift the right level e T = — V /2 + eV g by changing 
the gate voltage V g at the right dot, similar to the experiment |4j. Due to capacitive 
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Table 1. Parameter values used for the numerical calculations, extracted from the 
experiment Hj. 



parameter 


value 


capacitive level offset 


Se 


= 30/ieV 


interdot tunnel-coupling 


T c 


= 4.5 fjjeV 


tunnel rates 


r L = r R 


= 30/ieV 


source-drain voltages 


V 


= 650 /iV 




V x 


= 0//V 


electron temperature 


1/0 


= 10/ieV 



couplings there is an additional offset Se during the pulse and the energy of the right dot 
is e r = — Vo/2 + eV g + Se. Thus, the level offset e = e r — £i is e = eV g and E\ = eV g + Se. 
In principle it is sufficient to specify either e or ei, since this fixes automatically the 
other one. However, for convenience we will use both notations in the following. 

3.1. Numerical Solution of the TLQME 



We have investigated the pump-probe scheme presented in Sec. |1.1| using the TLQME 
given by Eqs. (fS]) and The operators are represented in the basis {|0) , \i) , |r)}. The 
resulting system of differential equations is propagated by means of an auxiliary-mode 



expansion described in Appendix A In order to calculate the number of pulse-induced 



tunneling electrons the numerically determined current [Eq. (12)] was integrated over 
time and the stationary contribution was subtracted. Therefore, the following integral 
has to be calculated 

+oo 

N(e, r) = - [ dt [J R (e, t) - J (e)] + . (16) 



The dependence on e accounts for both situations e = Eq and e = e\, which have a fixed 
relation E\ = e + Se as explained above. J refers to the stationary current reached 
at the end of the initialization and the measurement phase. The last term has been 
added to correct for the different stationary current during the manipulation phase (cf. 
Fig#). 

Figure |2^l shows the current Jo at the end of the initialization phase as a function 
of the level offset sq. This current is maximal when the energies of the two charge 
states coincide (eq = 0). Also shown is the analytical result for large source-drain 



voltages 34 35]. In Fig.[2p the number of pulse-induced electrons is shown as a 
function of pulse length r and energy difference e /i- Hereby, the function N(ei,t) is 
asymmetric around e\ = 0. Moreover, one can clearly observe an oscillatory behavior 
of N . The frequency of the oscillations increases with increasing values of e\. This 
is explicitly shown in Fig.[2b, where the function N(e, r) is shown for two energy 
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e 1 bieV] 



Figure 2. Numerical results for J and N as functions of energy difference e /i 
and pulse length r. a) Stationary current Jo(eo) before pulse. Symbols indicate 
numerical expressions and dashed lines denote analytical expressions for the LBL, i. e. 
Vsd °o 34 35] . b) Number of pulse-induced tunneling electrons N calculated with 
Eq. (16). The dashed white lines show expected positions of maxima r n = (2n+l)7r/f2 
with n — 0, 1, . . . c) Number of pulse-induced tunneling electrons N vs pulse length 
for constant energy difference Eq — und E\ = 0. 



differences, 6$ = and E\ = 0, respectively. For the latter case, the oscillations have the 
largest amplitude. 

Following the hypothesis that the number of tunneling electrons reflects just the 
occupation of the right dot at the end of the voltage pulse, one may naturally interpret 
the oscillations seen in Fig.[2|D as Rabi oscillations between the two charge states [2], 
with the frequency given by 

n= vV + 4T c 2. (17) 

Consequently, the maxima of iV should appear at pulse lengths r n = (2n+l)ir/Q with 
n = 0,1,..., whereby the frequency depends on E\ according to Eq. (fl7|). These 
positions are indicated by white dashed lines in Fig.[2]D. Around E\ = the maxima of 
iV are indeed found at the expected positions. However, for large E\ the maxima are 
clearly shifted compared to r n . 

It is interesting to compare the numerical data for N(e, t) with the instantaneous 
occupation of the right dot o~ IT (t). This is shown in Fig.[3] for two energies e. For 
£i = (e = —30 fieV) both quantities are almost identical to each other. In the 
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initialization 



manipulation 



measurement 




time, pulse -length [ns] 

Figure 3. Occupation of the right quantum dot cr rr (red lines) and the number of pulse- 
induced electrons N (blue lines) as functions of time t or pulse length t, respectively, 
from solving the TLQME numerically. Maxima are expected at t n — r n = (2n+l)7r/fi 
with n = 0, 1, . . . (vertical dashed lines). 



case £i = 30/ieV (eo — 0) one observes considerable deviations between N and a„. 
In contrast to the transferred electrons N, the occupation cr rr does have maxima at 
the times expected from the Rabi oscillations. Furthermore, N also takes negative 
values. Altogether, this is in clear contradiction to the assumption that the number of 
transferred electrons is just the occupation of the right dot. 

In summary, the numerical results and the experimental observations are 
qualitatively very similar. In particular, the described features of the behavior of N 
can also be seen in the experimental results. However, the reasons for the peculiar 
features, the asymmetry of iV around e\ = 0, the negative values of N for large E\ and 
the apparent shift of the maxima, cannot be explained by the numerical investigations 
alone. Therefore, in the remaining part of this section we will consider an analytically 



solvable model based on the Markovian quantum-master equations of Sec. 2.3 



3.2. Model Calculations in the LBL 



Starting point of the analysis is the initialization phase described by Eq. (13). Its 



stationary state provides the input for the manipulation phase. The stationary solution 



of Eq. (13) directly gives the the initial state within the considered pump-probe scheme. 
Setting a = yields the stationary populations au, a TT and coherences o& 



Mo) 

(7 rr (0) 



1 - 



8T 2 



12T C 2 + T 2 + 4e 2 
4T 2 



2 - 



12T C 2 + T 2 + 45 2 ' 



(18a) 
(18b) 
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2iT c (r + 2i g0 ) 

^ r(0) = " + P + • (18C) 

The stationary state is in general a mixed state. Considering the case of strong coupling 
to the reservoirs, r>T c , gives 

a a {0) = 1 - 0{T 2 C ) , (19a) 

MO) = 0{Tj) , (19b) 
2iT 

MO) = - = J- + 0(T C 2 ) . (19c) 

t R — ^1^0 

Obviously, in this case the initialization leads to an almost perfect localization of an 
electron in the left quantum dot. However, it is important to notice that at the same time 
finite coherences are unavoidable, which will have consequences for the manipulation 
phase. 

For the manipulation phase it is convenient to introduce the following combinations 
of matrix elements of the density matrix a, 

s(t) = (T U + a„ , w(t) = (T U - (T rr , 

u(t) = (x & + cx rf , v(t) = -i (<Jfy — Pre) ■ (20) 



The respective equations of motion [Eqs. (14)] can be solved using a Laplace 



transformation 36 , which is done in Appendix B Under the assumption i\ = I\ 



one finds for the total occupation 

+ 7 + 2r 



9F 

s(t) = e-^ +ar )*s(0) + [l- e -(T +2r ) t ] . (21a) 



and for the difference 

w (t) = ^$ 1 e-~< t [1 - cos(flt)] u(0) 

+ ^e^* sin(fit)^(0) + £ i +4T c ^M e -%(0) . (21b) 

Together, these two expressions allow for a calculation of the occupations in the DQD. 
Moreover, the current through the right barrier is given by 

J for t < , 

J R (t)={ e^a„(t) -eV[l-s(t)} for 0<t<r, (22) 

eT a VI (t) for r < t . 

Hereby, the stationary current is J = ercr rr (0) = er<r rr (oo). The typical time 
dependence of Jr(£) for 7 = is shown in Fig.jlJ. 

Thus, we are ready to get the number of pulse-induced tunneling electrons 



according to Eq. (16). This total number has contributions from the manipulation and 



measurement phase, i. e. N(e, r) = Ni(e, r) + N (e, r), which are according to Eq. (22) 



Ni(e,r)= [ dt{ 7 t7 ir (*)-r[l-s(t)]} (23a) 
Jo 

dt [r<T a (t) - Jo/e] . (23b) 
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The first contribution can explicitly be obtained from Eq. (21). For the case 7 = and 
s(0) ~ 1, one finds s{t) « 1 and consequently 

N 1 (e,r)^0. (24) 

For the second contribution we use again a Laplace transform, cf. Appendix B| The 
result of this procedure is the following expression for the number of pulse-induced 
electrons in the measurement phase, 

N Q (e, t) = X - [u 2 (0) - v 2 (0)] - 20-«(O)0 ir (O) (25) 

- - [u(0)u(r) - v(0)v(t)} + a a {0)a a (r) + a«(0)a rr (r) . 

Obviously, the value of A"o depends in a non-trivial way on all elements of the reduced 
density matrix. For weak inter-dot couplings, T c <C T, we can use Eq. (19) for the initial 
density matrix. As in Eq. (19) we keep only terms linear in T c /Y. Thereby we obtain, 
together with Eq. (24), a compact expression for the number of pulse-induced electrons 

(26) 



N(s,t) 



(r) 



AT r 



T 2 + Ael 



£qu(t) + 2 v(t) 



This is a central result of this work as it shows quantitatively the differences between 
A" and a rx and, as we discuss in the following, explains the main features seen in the 
experiment. 

Firstly, Eq. (26) implies that for large energy differences, e§ 3> T c , the number of 
pulse-induced tunneling electrons indeed corresponds to the occupation of the right dot 
at the end of the pulse. Therefore, it is a good measure for the occupation only for 
non-resonant initialization as already seen in Fig.[3j In this case the second term in 
Eq. (26) becomes very small. However, if the energy difference Eq, which is relevant in 
the initialization and in the measurement phase, vanishes, this term cannot be neglected. 
The corrections may lead to A" taking negative values, which is shown in Fig.|4j Secondly, 
for energy differences E\ close to zero one finds N(e,t) ~ er rr (r) = |[s(r) — w(r)] ~ 
\[1 — w(t)]. In this case, A" is mainly determined by the population difference w(r), 
which is explicitly given by Eq. (21b). An important consequence of this dependence is 
the occurrence of an asymmetric behavior of w(t) as a function of e± for a finite real 
part u(0) of the initial coherences. For example, considering times t n = (2n+l)7r/f2 
with n = 0, 1, . . . and assuming for the moment 7 = 0, Eq. (21b) yields 



w (t n ) 



w(0) 2 47>(0) 4T c 2 w(0) 



£1 - 



(27) 



n 2 1 n 2 n 2 

Obviously, this expression is only symmetric with respect to the energy difference E\ 
if u(0) = 0. In general, one finds w(t n ;£i) 7^ w (t n ; — Ex). In summary, the analysis 
within the described model provides a good explanation for the main features seen in 
the numerical results. These features result from an additional dependence of A^ on all 
the matrix elements of the reduced density matrix. 

Finally, we will briefly discuss damping of the oscillations. To this end, the 
analytical result given by Eq. (23) and the numerically obtained results from Fig. [3] are 
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Figure 4. Number of pulse-induced electrons N as a function of pulse length r. 
The numerically obtained results (squares) from Fig. [3] are shown together with the 
analytical expression [Eq. (23)] of the LBL model for two damping coefficients 7 (lines). 



shown together in Fig.El Without additional damping processes, i. e. 7 = in Eqs. (14) 



one finds for both energy differences undamped oscillations (dashed lines in Fig. [4]). 
Hereby, the frequency and the positions of the maxima are in good agreement with 
the numerical results. For e\ = 30 /xeV the number of pulse-induced electrons, N, takes 
negative values. An almost perfect agreement with the numerical results can be achieved 
by introducing decoherence into the manipulation phase (7 > 0). Taking 7 = 0.46 /xeV 
(full lines in Fig. [5]) yields a very good description of the pulse- length dependence of N. 
Obviously, decoherence is a necessary ingredient to obtain a consistent picture. In the 
numerical calculations, the decoherence originates from the finite source-drain voltage 
and the broadening of the energy levels. This leads to a non-vanishing probability of 



the electron leaving the DQD 37 . In the experiment other sources of decoherence exist 



and these will typically dominate the damping of coherent effects. The most important 
processes in this regard are interaction of electrons with phonons, background-charge 
fluctuations and cotunneling jlj. From the energy-difference dependence of the damping 
rate, which is extracted from N(s,t), one may infer about the nature of the relevant 
decoherence processes (3l|4j. 
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4. Conclusions 

In summary, we have theoretically investigated a pump-probe scheme realized in a 
recent experiment on the coherent manipulation of charge qubits in double quantum 
dots [I]. To this end we have numerically simulated the pump-probe scheme using a 
time-local quantum master equation. The equations are solved by the auxiliary-mode 
expansion technique described in |Appendix A which, in general, provides a flexible 



and viable method to study time-resolved electron transport. The numerical results for 
the number of pulse-induced electrons N(e, r) show good qualitative agreement with 
the experimental results. In particular, the main feature seen in the experiment, i.e. 
clear oscillations of iV as a function of pulse length r, are also observed. Moreover, 
two other characteristics of the experimental results are seen in the numerical data: 
the asymmetry of A" around E\ — and the occurrence of negative values. To adress 
these, so far unexplained, features we considered the Markovian limit of the QME in 
the respective pump-probe stages. The resulting equations can be solved analytically 
and lead to the main result of the article, namely the expression for the transferred 
electrons N. It turns out, that the value of A" depends in a non-trivial way on all 
elements of the reduced density matrix. Only for large initial energy offsets, £o > T c , 
the number of pulse-induced tunneling electrons corresponds to the occupation of the 
right dot at the end of the pulse. For small values of Eq, larger differences between 
the occupation and the number of pulse-induced electrons are expected, which explains 
the occurrence of negative values. Further analysis shows that the asymmetry can be 
attributed to unavoidable initial coherences resulting from the initialization stage. Thus, 
the presented findings strengthen the common interpretation of the observed features 
in terms of coherent manipulation. 

The expression for N(e, t) may readily be tested against the experimental results. 
An interesting question for further studies would be the analysis of the decoherence rate 
as a function of energy difference. To this end, one could extent the numerical description 



by including phonons 38 , 39 or background-charge fluctuations. The decoherence rate 
may then be extracted from the numerical data using the analytical expression for the 
pulse-induced tunneling electrons from the LBL model. 

Appendix A. Auxiliary Mode Propagation 



In order to perform the energy integration in Eqs. (10) we expand the Fermi function 
/(e) as a finite sum over simple poles 

with xt = H±x P /(3 and Imx p > 0. Here, [i is the chemical potential and ft the 



inverse electron temperature. Instead of using the Matsubara expansion 40 , with poles 
x p = br(2p— 1), we use a partial fraction decomposition of the Fermi function ETl, which 
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converges much faster than the standard Matsubara expansion. In this case the poles 

The poles are 



x 



±2^/1^ are given by the eigenvalues z p of an N^xNp matrix 



41 



arranged such that all poles xt (x P ) are m t ne upper (lower) complex plane. 



Employing the expansion given by Eq. (A.l), one can evaluate the energy integrals 



necessary to compute the reservoir correlation functions given by Eqs. (10) by contour 



integration. Thereby, the integrals in Eqs. (10) become (finite) sums of the residues. 
One gets for t > t' 



N F 
p=l 



>/ t U"xS P (t") 



with the auxiliary modes for reservoir a given by 

Xap{t) = \p-a + A Q (t)] ±x p /(3± iT amm /2. 



(A.2a) 
(A.2b) 

(A.2c) 



Hereby A a (t) is due to the time-dependent single-particle energies e a k(t) of the reservoir 
Hamiltonian [Eq. (pjj)]. The last term containing T amm leads to the correct broadening 



of the energy levels as discussed in Sec. 2.2 



With the expansion (A. 2) one can write the the auxiliary operators [Eq. ([9 
terms of the partial correlation functions. One readily obtains 

1 



in 



I - q(x) 

L amm u rn 



(A.3) 



with the "partial" auxiliary operators 



(A.4) 



to 



An analogous expression holds for Aam- The partial auxiliary operators A«m ;p contain 



(i) 



instead of the full correlation function just an exponential factor (see Eq. (A. 2b)) and 
the corresponding equation of motion, 

J^SU*) = l f«m m s£ x) -i [#s(*),a£UC 

+ i*xL(*)A£> (t) , (A.5) 



is easily propagated with the initial values Aam;p{to) = 0. 

For the numerical calculation in Sec. [3] we represent all operators in the basis 



{|0) , \£) , |r)}. This renders the equations of motion, Eq. ([8]) and Eqs. (A.5), to become 
matrix equations. However, note that these ~ 2A^p equations of motion (we use A^p=80) 



are not coupled. Moreover, the density matrix o does not enter Eqs. (A.5), which makes 



the propagation scheme very efficient. In order to propagate the equations we use a 
fourth order Runge-Kutta scheme 42 with a constant time-step 8t = 0.0002 //zeV. 
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Appendix B. Analytical Solution of the LBL Model 

We summarize the solutions of the density matrix equations given in Sec. |2.3| and used in 



Sec. 3.2 They are most easily obtained using Laplace transforms x(s) = J °° dte st x(t). 



Manipulation Phase 

Assuming T L = T R = T renders the equation for the total occupation s independent 



from the others, s(t) = — [j+2T]s(t) + 2T, with the solution given by Eq. (21a). The 
other components w(t), u(t), and v(t) are given by coupled differential equations which 
are solved by applying a Laplace transform. The resulting linear system reads 

/ s+7 -2 



V 



o 

2T r 



s+7 
-ei 





( w(s) 


\ 


/ 


w(0) 


\ 




u(s) 






«(0) 




) 


{ w(s) 


/ 


V 


v(0) 


J 



(B.l) 



By means of the solutions w(s), u(s), and v(s) one gets with an inverse Laplace transform 
the final expressions 

- 4T c 2 cos(^) +g ^ (o) + 2r c£l [l-cos(m)] + 2T c sin(nt) - 

)2 Q2 (\ 



W(t) 
U (t) 

v(t) 



- 7 t 



fi 5 



n 



2T c e x [l-cos(fit)] + 4r 2 +e 2 cos(fit) ^ (o) _ gisin(Ot) 



fi 2 

2T C sin(fit) 

fi 



w(0) 



£i sin (fit) 
fi 



fi 2 

u(0) + cos(fit)w(0) 



fi 



(B.2a) 
(B.2b) 
(B.2c) 



where we have used the Rabi frequency fi from Eq. (17). 

Measurement Phase 



The Laplace transformation of Eqs. (13) leads to the following linear system 
/ s+T r 



V 





-iT c 
iT c 



s+T 
-iTc. 



-iT c 
iT c 
s — ie H- ' 




-iT c 




\ 



s+ie +2 / 



/ o u {s) \ 

CT rr (s) 

a ir (s) 



\ 



cr TT {T) 

<7£r(r) 



(B.3) 



This system may be solved for the components au(s) , <7 rr (s) , a£ T (s) and a T i(s). In order 



to calculate the number of pulse-induced tunneling electrons [Eq. (25)], one just needs 



cr rr (s). The stationary value of cr rr for t — > oo can be calculated through the following 



limit, 



lim s cx rr (s) 



4T 2 



12T 2 +r 2 +4e 2 



<Jrr OO 



(B.4) 



This expression is identical to Eq. (18). The integrated current is given by 

4 



limT [cr rr (s) — cr rr (oo)/s] 

s— >0 



T 2 c a u {r) + [T 2 +(r/2) 2 +£ 2 ] a rr ( 



12T 2 +r 2 +4e 2 

- 2T C [e Re a ir (r) + (T/2) Im a h {r)\ 

- 4T 2 [2T 2 +r 2 ] / [l2T 2 +r 2 +4 £ 2 ] 



(B.5) 
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Using the expressions for the stationary density matrix [Eqs. (18)] finally yields Eq. (25). 
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